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ABSTRACT 

Spectral lag, which is defined as the difference in time of arrival of high and low energy pho- 
tons, is a common feature in Gamma-ray Bursts (GRBs). Previous investigations have shown 
a correlation between this lag and the isotropic peak luminosity for long duration bursts. How- 
ever, most of the previous investigations used lags extracted in the observer-frame only. In this 
work (based on a sample of 43 Swift long GRBs with known redshifts), we present an anal- 
ysis of the lag-luminosity relation in the GRB source-frame. Our analysis indicates a higher 
degree of correlation —0.82 ± 0.05 (chance probability of ~ 5.5 x 10~ 5 ) between the spectral 
lag and the isotropic peak luminosity, L lso , with a best-fit power-law index of — 1.2 ± 0.2, 
such that Li SO oc lag -1 ' 2 . In addition, there is an anti-correlation between the source-frame 
spectral lag and the source-frame peak energy of the burst spectrum, E p ^(l + z). 

Key words: Gamma-ray bursts 



1 INTRODUCTION 

Gamma-ray Bursts (GRBs) are extremely energetic events and pro- 
duce highly diverse light curves. A number of empirical correla- 
tions between various properties of the light curves and GRB en- 
ergetics have been discovered. However, the underlying physics of 
these correlations is far from being understood. 

One such correlation is the relation betwe en isotropic peak lu - 
minosity of long bursts and their spectral lags dNorris et al. Il200(jh . 
Various authors have studied this relation us ing arbitrary observer- 
frame energy-bands of various instruments jUkwattaetalJl2010j; 
Hakkila et all 120081 ; ISchaefeJ 120071 : iGehrels et alj 120061 ; iNorrisl 
2002}) . These investigations support the existence of the relation, 
howev er with considerable scatter in the extracted results. Re- 
cently, iMargutti et al.1 fcOld) investigated spectral lags of X-ray 
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flares and found that X-ray flares of long GRBs also exhibit the 
lag-luminosity correlation observed in the prompt emission. 

The spectral lag is defined as the difference in time of arrival 
of high and low energy photons and is considered to be positive 
when the high-energy photons arrive earlier than the low energy 
ones. Typically the spectral lag is extracted between two arbitrary 
energy bands in the observer frame. However, because of the red- 
shift dependence of GRBs, these two energy bands can correspond 
to a different pair of energy bands in the GRB source-frame, thus 
potentially introducing an arbitrary energy dependence to the ex- 
tracted spectral lag. 

In order to explore whether the lag-luminosity relation is in- 
trinsic to the GRB, it is preferable to extract spectral lags in the 
source frame as opposed to the observer frame. At least two cor- 
rections are needed to accomplish this: 1) Correct for the time dila- 
tion effect (z-correction), and 2) Take into account the fact that for 
GRBs with various redshifts, observed energy bands correspond 
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to different energy b ands at the GRB source-frame (K-correction; 
iGehrels et alj fcood) ). 

The first correction is straightforward and is achieved by mul- 
tiplying the extracted lag value (in the observer- frame) by (1 + 
z)^ 1 . T he second correction, on the other hand, is not so straight- 
forward. IGehrels et all d2006l) attempted to approximately correct 
the spectral lag by multiplying the lag value (in the observer- 
frame) by (1 + z) ' 33 . We note here that this correction is based 
on the assumption that the spectral lag is proportional to the pulse 
width and that the pulse width itself is pr oportional to the energy 
dZhangetalJl2009l ; iFenimore et aTlll995h . These approximations 
depend on clearly identifying corresponding pulses in the light 
curves of each energy band, and may be of limited validity for a 
large fraction of GRBs in which the light curves are dominated by 
overlapping multi-pulse structures. 

Using a sample of 31 Swift GRBs. lUkwatta et alj falOal) 
(hereafter U10) found that the correlation coefficient improves sig- 
nificantly after the z-correction is applied. However, this correlation 
does not imp rove further after the application of the K-correction 
as defined by IGehrels etal] J2006h ■ 

An alternative is to make the K-correction by choosing two 
appropriate energy bands fixed in the GRB source-frame and pro- 
jecting these bands into the o bserver-frame using th e relation 

server = E s 

ourcc 

/(l + z). lUkwatta et all feOlObl) used this 
method for the first time to investigate the lag-luminosity relation 
in the source-frame of the GRB. They selected two source-frame 
energy bands (100 - 200 keV and 300 - 400 keV) and used back- 
ground subtracted as well as non-background subtracted Swift 
data to extract lags. Non-background subtracted data were used to 
improve the signal-to-noise ratio for weak bursts. They found that 
the source-frame relation se ems a bit tighte r , but w ith a slope con- 
sistent with previous studies. lArimoto et al.l d20l5T) also looked at a 
limited sample of HETE-II bursts (8 GRBs) both in the observer- 
frame and the source-frame and concluded that there is no signif- 
icant effect from the redshift. However, the redshift distribution of 
th eir burst sample is very narrow and peaks around one. In contrast 
to lUkwatta et alj fcOlObh . in this study we used only background 
subtracted data and measured the lag between source-frame energy 
bands 100 - 150 keV and 200 - 250 keV (the reason for selecting 
these particular energy bands is described in §[2j( for a sample of 43 
Swift bursts with spectroscopic redshifts. 

In this work we have investigated only long GRBs, i.e., bursts 
with duration greater than ~ 2 seconds. It is rather difficult to 
test the Lag-Li so relation effectively for short GRBs due to a lack 
of spectroscopically measured redshifts. None of the short bursts 
detected so far have any redshift measurements obtained from 
a spectroscopic analysis of their optical afterglow. Moreover, it 
has been shown that short GRBs have either s mall or negligible 
lags dNorris & Bonnellll2006l : IZhang et alj|2006h - According to the 
Lag-Li so relation, these small lag values imply that short bursts to 
be highly luminous. However, based on the redshift measurements 
of their host galaxies we can show that short GRBs are generally 
less luminous than long bursts. Hence short bursts se em to not fol- 
low the lag-luminosity relation dGehrels et ailj2006l) ). 

The structure of this paper is the following: In § [2] we dis- 
cuss briefly our methodology for extracting spectral lags. In §[3]we 
present our results for a sample of 43 Swift GRBs. We discuss our 
results with two candidate models in §|4] Finally, in the last section 
(§ O we summarize our results and conclusions. Throughout this 
paper, the quoted uncertainties are at the 68% confidence level. 



2 METHODOLOGY 

The Swift Burst Alert Telescope (BA T) is a highly sensitiv e 
instrument using a coded-mask aperture dBarthelmy et al.|[2005h . 
BAT uses the shadow pattern resulting from the coded mask to fa- 
cilitate localization of the source. When a gamma-ray source illu- 
minates the coded mask, it casts a shadow onto a position-sensitive 
detector. The shadow cast depends on the position of the gamma- 
ray source on the sky. If one knows the tile pattern in the coded 
mask and the geometry of the detector, it is possible to calculate 
the shadow patterns created by all possible points in the sky using a 
ray-tracing algorithm. Hence, by correlating the observed shadow, 
with the pre-calculated shadow, one can find the location of the 
source. However, each detector can be illuminated by many sources 
and a given source can illuminate many detectors. Hence, in or- 
der to disentangle each sky position, special algorithms have been 
developed and integrated in to the data analysis software by the 
Swift BAT team. 

To generate background-subtracted light curves we used a pro- 
cess called mask weighting. The mask weighting assigns a ray- 
traced shadow value for each individual event, which then en- 
ables the user to calculate light curves or spectra. We used the 
batmaskwtevt and batbinevt tasks in FTOOLS to gener- 
ate mask weighted, background-subtracted light curves, for various 
observer-frame energy bands, as shown in Table Q] These are the 
energy bands that correspond to fixed energy bands in the source- 
frame i.e. 100 — 150 and 200 — 250 keV. These particular en- 
ergy bands were selected so that after transforming to the observer- 
frame they lie in the detectable energy range of the Swift BAT 
instrument (see Fig. [TJ. Even though the BAT can detect photons 
up to 350 keV, we limited the upper-boundary to 200 keV in the 
observer-frame. This is because the mask-weighted effective area 
of the detector falls rapidly after 200 keV and as a result the contri- 
bution to the light curve from energies greater than ~200 keV (in 
observer-frame) is negligible dSakamoto et al.l201 lh . 

The spectral lags were extracted using the improved cross- 
correlation function (CCF) analysis method described in U10. In 
this method, the spectral lag is defined as the time delay corre- 
sponding to the global maximum of the cross-correlation function. 
The CCF with a delay index d is defined as, 

Emin(JV,JV-d) 
i=max(i,i-d) a:i 2'*+ d 
UUK(d,x,y) = - (1) 

v Ei x i E, Vi 

where Xi and yt are two sets of time-sequenced data spread over 
N bins. The time delay is obtained by multiplying d by the time 
bin size of the light curves. A Gaussian curve was fitted to the CCF 
(plotted as a function of time delay) to extract the spectral lag. The 
uncertainty in the spectral lag is obtained by simulating 1,000 light 
curves using the Monte Carlo technique (see U10 for more details). 

The isotropic peak luminosity (L lso ) and its uncertainty for 
each GRB is obtained using the method described in U10. In 
essence, a typical GRB sp ectrum can be described by the Band 
function dBandetal.lll993l) , for the photon flux per unit photon en- 
ergy using 

N(E) _ J M^re-^^, E * (^)E pk 

[ ^HlOOkoW l(2+a)100keVJ c , cioc, 

which has four model parameters: the amplitude (A), the low- 
energy spectral index (a), the high-energy spectral index (J3) and 
the peak (Ep^) of E 2 N(E) spectrum (also called the vF v spec- 
trum, apart from a factor of Planck's constant). Using these spectral 
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Figure 1. Fixed energy bands at the GRB source-frame are projected to various energy bands at the observer-frame, depending on the redshift. 



parameters the observed peak flux can be calculated for the source- 
frame energy-range Ei = 1.0 keV to E2 = 10, 000 keV using 

/•S 2 /(l+*) 

/obs= / N{E)EdE. (3) 

Je 1 /(1 + z) 

The isotropic peak luminosity is defined by 

L iso = Audi f obs (4) 
where dh is the luminosity distance: 

dL =^r , d " ( 5 ) 

For the current universe we take Qm = 0.27, Q. l = 0-73 and 
the H ubble constant Ho to be 70 (kms _1 )/Mpc dKomatsu et al.1 
l2009h . For more details of the L; so calculation see U10. 



3 RESULTS 

We employed an additional 12 long bursts to the GRB sample (31 
GRBs) that was used in U10, which increased the total sample to 
43. This sample has redshifts ranging from 0.346 (GRB 061021) to 
5.464 (GRB 060927) with an average redshift of ~2.0. The spectral 
information for the additional 12 bursts used in this paper is given in 
Table|2] The calculated peak isotropic luminosities, spanning three 
orders of magnitude, are given in U10 and Table|2] 

By choosing appropriate energy bands in the observer- 
frame (according to the redshift of each burst), we extracted 
mask-weighted background-subtracted light curves for the se- 
lected source-frame energy bands 100-150 and 200-250 keV. The 
observer-frame energy bands used for each burst are shown in Ta- 
ble Q] Note that the energy gap between the mid-point of the two 
source-frame energy bands is fixed at 100 keV whereas in the 
observer-frame, as expected, this gap varies depending on the red- 
shift of each burst (see the TableQJ. For example, in GRB 060927, 
this gap is 16 keV and in GRB 061021, it is 75 keV. This is in con- 
trast to the spectral lag extractions performed in the observer-frame 
where this gap is treated as a constant. 

The extracted spectral lags for the source-frame energy bands 



100-150 and 200-250 keV are listed in Tabled The Swift BAT 
trigger ID, the segment of the light curve used for the lag extraction 
(T + Xs and T + Xe, T is the trigger time), the time binning of 
the light curve, and the Gaussian curve fitting range of the CCF 
vs time delay plot (with start time, and end time denoted as LS 
and LE respectively) are also given in Table[5] Of 43 bursts in the 
sample there are 24 bursts which have lags greater than zero. The 
remaining 19 bursts either have lags consistent with zero (16 bursts) 
or negative values (3 bursts). 

For the 24 bursts which have positive lags with significance 
1-cr or greater (see Table[3), we find that the redshift corrected lag 
is anti-correlated with Liso- The correlation coefficient for this re- 
lation is -0.82 ± 0.05 with a chance probability of ~ 5.54 x 10~ 5 . 
The extracted correlation coefficient is significantly higher than the 
correlation coefficient (averaged over the six combinations of stan- 
dard BAT energy channels) of ~ —0.68 reported in U10. Various 
correlation coefficients of the relation are shown in Table [4] where 
uncertainties in the correlation coefficients were obtained through 
a Monte Carlo simulation utilizing uncertainties in Li so and the lag 
values. The null probability that the correlation occurs due to ran- 
dom chance is also given for each coefficient type. 

Fig. [2] shows a log-log plot of isotropic peak luminosity vs 
redshift-corrected spectral lag. The solid line shows the following 
best-fit power-law curve: 



log (^-) = (54.7 ± 0.4) - (1.2 ± 0.2) log La 9^\ (6 ) 
\ erg/s J 1 + 2 



Since there is considerable scatter, the uncerta inties of the 
fit parameters are multiplied by a factor of -J x 2 /ndf = 

\J 84.36/22 = 1.96. The dash lines indicate the estimated 1-cr 
confidence level, which is obtained from the cumulative fraction of 
the residual distribution taken from 16% to 84%. 

The best-fit power-law index (—1.2 ± 0.2) is con sistent with 

observer-frame results obtained bv lNorris et al.ld2000h ( 1.14) 

and the average power-law index of —1.4 ± 0.3 reported in U10. 
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Figure 3. All combinations of fixed observer-frame energy channel (canonical BAT energy bands: channel 1 (15-25 keV), 2 (25-50 keV), 3 (50-100 keV) 
and 4 (100-200 keV)) spectral lag values as a function of fixed source-frame energy channel (between 100-150 keV and 200-250 keV) lag values. Black and 
red data points and labels corresponds to redshift uncorrected and corrected cases respectively. The blue dashed line corresponds to the equality line of the two 
parameters in each panel. 
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Table 1. The observer-frame energy bands and energy gaps (the energy difference between the mid-points 
of energy bands) for bursts in the sample. 



GRB 


Redshift 


Low Energy Band (keV) 


High Energy Band (keV) 


Energy Gap (keV) 


GRB050401 


2.899 1 


26-38 


51-64 


26 


GRB050603 


2.821 2 


26-39 


52-65 


26 


GRB050922C 


2.199 3 


31-47 


63-78 


32 


GRB051111 


1.549 4 


39-59 


78-98 


39 


GRB060206 


4.056 5 


20-30 


40-49 


20 


GRB060210 


3.913 6 


20-31 


41-51 


21 


GRB060418 


1.490 7 


40-60 


80-100 


40 


GRB060904B 


0.703 8 


59-88 


117-147 


59 


GRB060908 


1.884 9 


35-52 


69-87 


35 


GRB060927 


5.464 10 


15-23 


31-39 


16 


GRB061007 


1.262 11 


44-66 


88-111 


45 


GRB061021 


0.346 12 


74-111 


149-186 


75 


GRB061121 


1.315 13 


43-65 


86-108 


43 


GRB070306 


1.496 14 


40-60 


80-100 


40 


GRB071010B 


0.947 15 


51-77 


103-128 


52 


GRB071020 


2.145 16 


32-48 


64-79 


32 


GRB080319B 


0.937 17 


52-77 


103-129 


52 


GRB080319C 


1.949 18 


34-51 


68-85 


34 


GRB080411 


1.030 19 


49-74 


99-123 


50 


GRB080413A 


2.433 20 


29-44 


58-73 


29 


GRB080413B 


1.101 21 


48-71 


95-119 


48 


GRB080430 


0.767 22 


57-85 


113-141 


56 


GRB080603B 


2.689 23 


27-41 


54-68 


27 


GRB080605 


1.640 24 


38-57 


76-95 


38 


GRB080607 


3.036 25 


25-37 


50-62 


25 


GRB080721 


2.591 26 


28-42 


56-70 


28 


GRB080916A 


0.689 27 


59-89 


118-148 


59 


GRB081222 


2.770 28 


27-40 


53-66 


26 


GRB090424 


0.544 29 


65-97 


130-162 


65 


GRB090618 


0.540 30 


65-97 


130-162 


65 


LrKB090/ 1 JD 


O.UOU 


25-38 


50-63 


25 


GRB090812 


2.452 32 


29-43 


58-72 


29 


GRB090926B 


1.240 33 


45-67 


89-112 


45 


GRB091018 


0.971 34 


51-76 


101-127 


51 


GRB091020 


1.710 35 


37-55 


74-92 


37 


GRB091024 


1.091 36 


48-72 


96-120 


48 


GRB091029 


2.752 37 


27-40 


53-67 


27 


GRB091208B 


1.063 38 


48-73 


97-121 


49 


GRB100621A 


0.542 39 


65-97 


130-162 


65 


GRB 1008 14A 


1.440 40 


41-61 


82-102 


41 


GRB100816A 


0.800 41 


56-83 


111-139 


56 


GRB 100906 A 


1.727 42 


37-55 


73-92 


37 


GRB 11 0205 A 


2.220 43 


31-47 


62-78 


31 


GRB110213A 


1.460 44 


41-61 


81-102 


41 


(1)1 Watson et al 


(2003); (2) 


Berger & Beckeil |2005bl); (3) IPiranomonte et al l 42005 


\): (4) IPenprase et al. 



(9)lFvnboetalJi 



. . jFvnbo et alj 12003) : (11) Fynbo et al. (2009); (12) Fynbo et al. (2009J; (13) 
Fynbo et al J J2009h : (14) Ijaunsen et al l fcOOSl) ; ( 1 5) ICenko etalTHo07l): C16) ljakobsson* et alj J2007I) : (17) 
D'Eli a et all J2009t) ; (18 ) iFvnbo et all i2009h; ( 19) iFynbo et alj 120091); (2oT lFvnbo et al l i7m% ; (21) 



L> Elia et alj J2U0y); (IS ) lEvnbo et alj <2ULft1); ( 19) Ihvnbo et alj 120U9I); (.20) fynbo et al. 
Fvnbo et al I<2009l): (2 2) ICucchiara & Foxl J2008t) : ( 23)lFvnbo et all 120091) : ( 24)lFvnbo et all 



Prochaska et alj 


2009' 


Chornock et alJ 


2009) 



b; (26) 
;(30) 



Fvnbo et al 



^,2005 ); (25) 

120091) ; (27) lFvnbo et alJl2009l); (28 ) ICucchiara et alJl2008l): (29) 
120091); (31 ) JWiersema et alj l2009j); (32)|de Ugarte Postigo et aT 



Cenko et al . 

20091); (33) IFvnbo et all 120091); (34) Ichen etall 120091); (35) IXu et all 120091); (36) ICucchiara etal 



2009); (37) Chornock et al. (2009); (38) Wiersema et al. (2009); (39) Milvang-Jensen et al. (2010); (40) 



O'Meara et "alj J201fj); (41)lTanvir etall J20ld); (42)ICenko et aU 1201 ll); (43)lMilne & Cenkol J201 ll). 



4 DISCUSSION 

4.1 Spectral Lags: Observer-frame versus Source-frame 

U10 extracted spectral lags in fixed energy bands in the observer- 
frame and in this work for the same sample of 31 bursts we ex- 



tracted lags in fixed energy bands in the source-frame. In the 
observer-frame case, there are four energy channels (canonical 
BAT energy bands: channel 1 (15-25 keV), 2 (25-50 keV), 3 (50- 
100 keV) and 4 (100-200 keV)), thus six lag extractions per burst. 
It is interesting to study to what degree these different lags cor- 
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Table 2. GRB redshift and spectral information. Note that uncertainties of parameters that are reported with 90% confidence level have been reduced to 
1-cr level for consistency. 



GRB 


Peak Flux a 




Q 


I3 C 


L iso erg/s 


Reference 



GRB090812 

GRB090926B 

GRB091018 

GRB091020 

GRB091024 

GRB091029 

GRB091208B 

GRB100621A 

GRB100814A 

GRB100906A 

GRB110205A 

GRB110213A 



3.60 ±0.13 
3.20 ±0.19 
10.30 ±0.25 
4.20 ±0.19 
2.00 ±0.19 
1.80 ±0.06 
15.20 ±0.63 
12.80 ±0.19 
2.50 ±0.13 
10.10 ±0.25 
3.60 ±0.13 
1.60 ±0.38 



17Z 
500± 
6lj 



°' z -156 

91±\ 

f 100 
-100 

no 

-10 

124t» 
1061^ 

onn+46 
■ 4 '"_46 

98± 4 



-1 03+ 04 
1 - U 'J_0.04 

-0 13+ - 04 

u - ±o -0.04 
, t-o+0.24 
1 - do -0.37 

-0 20 +0 ' 25 
u - zu -0.25 

-1 10+° 13 

-1 46+ ' 17 
1+0.04 
-0.04 
-,+0.08 
-0.08 
4+0.08 
-0.09 
-,+0.06 
-0.06 
,+0.09 
-0.09 
.1 44+0-03 
1 -^-0.03 



-1.443 
-1.70j 
-0.64 4 

-l.ioj 

-1.52" 1 



-,+0.16 
-0.16 
n +0.31 
-0.31 
-2.44+ , 15 , 

-0.15 
1-0.01 
-0.01 
,+0.31 
-0.31 
,+0.31 
-0.31 
,+0.29 
-0.12 
hi. 44 
-1.44 
,+0.08 
-0.06 
-,+0.19 
-0.13 
,+0.31 
-0.31 

-2 36+ ' 31 
Z -' 3D -0.31 



-2.50Z 
-2.36^ 



-1.70: 
-2.36j 
-2.36j 
-2.32 4 
-2.45^ 
-2.02 4 
-2.20j 
-2.36^ 



7 
5.22 
6.96 
2.81 



+ 1.95\ 
-0.87' 
+3.88 n 
0.82^ 
1.76\ 



o.r,*; 

+0.19n 
0.16' 
5 56+ 2 ' 43 ) 
+0.60 \ 
0.15-1 
+0.65\ 
0.09^ 
+0.83\ 
0.34 J 
1.13\ 
0.69' 
+1.23 \ 
-0.43 1 
+0.57^ 



1.67 
1.68 
2.55 
8.27 
4.90 

2 - 78 -0.20^ 
3.53_ n -n ) 



x 10 52 
x 10 51 
x 10 51 
x 10 52 
x 10 51 
x 10 52 
x 10 52 
x 10 51 
x 10 51 
x 10 52 
x 10 52 
x 10 51 



Pal' Shin et al. (2009); Baumgartner et al. (2009) 

Briggs (2009) : Baumgartner et al. (2009) 
Golenetskii et al. (2009) ; Markwardt et al. (2009) 

Chaplin (2009) ; Palmer et al. (2009) 
Golenetskii et al. (2009) ; Sakamoto et al. (2009) 
Barthelmv et al. (2009) 
McBreen (2009); Baum gartner et al. (2009) 
Golenetskii et al. (2010) ; Ukwatta et al. (2010) 

von Kienlin (2010) ; Krimm et al. (2010) 
Golenetskii et al. (2010) ; Barthelmv et al. (2010) 
Golenetskii et al. (2011) 
Foley (2011) ; Barthelmv et al. (2011) 



1-second peak photon flux measured in photons cm s 1 in the energy range 15 — 150 keV. 



b Peak energy, E p ^, is given in keV. 

Values in brackets indicates estimated high-energy photon 



ISakamoto et alj2009l) . 



index, /3, which is the mean value of the BATSE /3 distribution iKaneko et afll200d : 



Table 4. Correlation coefficients of the lag-luminosity relation. 



Coefficient Type 


Correlation Coefficient 


Null Probability 


Pearson's r 


- 0.82+0.05 


5.54 x 10~ 5 


Spearman's r B 


- 0.70±0.06 


1.49 x 10" 4 


Kendall's r 


- 0.50±0.05 


6.63 x 10~ 4 



relate with source-frame lags (between fixed source-frame energy 
channels 100-150 keV and 200-250 keV). In Fig. [3] we show all 
combinations of observer-frame lags as a function of source-frame 
lags. The red data points show lags with the time-dilation correction 
due to cosmological redshift and black data points show lags with- 
out the time-dilation correction. From Fig.[3]it is clear that all plots 
show some correlation both in the time-dilation corrected (shown 
in red) and time-dilation uncorrected (shown in black) cases. We 
note that the correlation coefficients are greater than 0.5 in time- 
dilation uncorrected cases where BAT channel 1 is involved in the 
lag extraction. In the time-dilation corrected case all plots show cor- 
relation coefficients greater than 0.5 except for the Lag43 plot. De- 
spite these moderate correlation coefficients, the large scatter seen 
in these plots indicate that the observer-frame lag does not directly 
represent the source-frame lag. 



4.2 Lag-Z/iso Relation: Observer-frame versus Source-frame 

There are two important changes in the Lag-Luminosity relation 
which may occur when going from fixed observer-frame energy 
bands to fixed source-frame energy bands: A change in the power- 
law index, and a change in the dispersion of the data measured by 
the correlation coefficient. Table [5] summarizes these two parame- 
ters for various energy bands both in the observer-frame and in the 
source-frame. 

In the observer-frame the power-law index varies from ~ 0.6 
to ~ 1.8, with mean around 1.3. In the source-frame the index 
changes from 0.9 to 1 .23 with a mean of ~ 1 . 1 . Meanwhile, the cor- 
relation coefficient varies from 0.60 to 0.79 in the observer-frame 
and in the source-frame it changes from 0.76 to 0.90. Hence, ac- 



cording to Table[5] the source-frame Lag-Z/i so relation seems to be 
tighter than the observer-frame case with a slope closer to one. 



4.3 Spectral Lag - Ep^ Relation 

Now we investigate the relation between source-frame spectral lag 
and source-frame average peak energy (E p k(l + 2)) of the burst 
spectrum. In Fig. [4] we plotted _E p k(l + 2) as a function of source- 
frame lags. There is a correlation between these two parameters 
with a correlation coefficient of —0.57 ± 0.14. Various correlation 
coefficients of the relation are shown in Table [6] with uncertainties 
and null probabilities. 

The best-fit is shown as a dashed line in Fig. [4] yielding the 
following relation between _E p k(l + 2) and Lag/(1 + z): 



log 



£pk(l + z) 

keV 



= (3.7±0.1)-(0.56±0.06) log 



Lag /(ms) 
l + z 



.(7) 



The uncertainties in the fitted parameters are expressed with 
the factor of ^x 2 /ndf = ^30.71/22 fa 1.18. 

According to equation (6), -Liso oc (Lag /(l + z))~ 1,2 . 
From t he Yonetoku relation we know that Li so oc (_E p k(l + 
z)) 2 l l Yonetoku et ai]|2004 . Hence, from these two relations we 
expect to see a correlation between E p ^(l + 2) and Lag/(1 + 2) 
such as i?pk(l + z) oc (Lag /(I + z))~ ' 6 . 

The best fit slope of 0.56±0.06 is consistent with the expected 
slope of ~ 0.6 based on the source-frame lag-luminosity and the 
Yonetoku relation. However, note that the correlation coefficient is 
significantly smaller than the coefficient for the lag-luminosity rela- 
tion. This lower degree of correlation may be suggestive of bright- 
ness and dete ctor related select ion effects that have been noted in 
the literature dButler et alj2007h for the Yonetoku relation. 



4.4 Some Models for Spectral Lags 

U10 and this work have provided more evidence for the existence of 
the lag-luminosity relation based on a sample of Swift BAT GRBs 
with measured spectroscopic redshifts. This analysis calls for a 
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Table 3. Source-frame spectral lag values of long duration Swift BAT GRBs 



GRB 


Trigger ID 


T + X$ (s) 


T + X E (s) 


Bin Size (ms) 


LS (s) 


LE (s) 


Lag Value (ms) 


Significanc 


VJIVDUJUt'U 1 


1 13120 


23 03 


29 43 


64 


-2 00 


2 00 


310±145 


2 14 


VJIYU UJV7UUJ 


1 3 1 SfiO 
i j i 


-3 83 


3 08 


16 


-0 40 


40 


-16±21 


-0 76 


V I J\ I J \J J\J~Z.Z.v^ 


156467 


-2 70 


2.94 


16 


-1.00 


1 00 


1 ^64-68 


2.00 


riRROSI 1 1 1 

VJIYUVJ J 1111 


163438 


-6.96 


28.62 


64 


-4.00 


4.00 


33 3 ±251 


1.33 


VJIVDUOUZVJU 


1 804 SS 


1 9Q 


8 18 


16 


-2 00 


2 00 


86il 1 1 


77 


nRR0609i o 


180977 


-3 37 


5 08 


128 


-4 00 


4 00 


6S84-9SQ 


2 54 


nRR06041 8 


^0585 1 


-7 66 


33 04 
J J . U4- 




9 00 


9 00 


1 1 n_|_ 1 06 


-1 04 


HRR060Q04R 

VJ 1\JJ \J U \J y Irt o 


998006 


-1 97 


10 32 


5 12 


-6 00 


6.00 


1 24±436 


28 


hrro^oqos 


998^81 


-10 91 


3 68 

J. Do 




9 00 


2 00 


78il24 


63 
U.U 3 


nRR060Q97 


231362 


1 6Q 


8 04 


39 

jZ 


1 00 


1 on 




24 


HRR06I 007 
vjiyuvju i i 


939683 


23 86 


65.08 


4 


-0.20 


20 


52±22 


2.36 


HRR06I 09 1 

VJ IYUVJU 1 UZ 1 


234905 


-0.46 


14 64 


512 


-4.00 


4.00 


-430±975 


-0.44 


(TRR061 191 


239899 


60 44 


80 66 


4 


-0 20 


20 


22±10 


2.20 


HRR070306 

VJIYUVJ / UJUU 


963361 

iUJJU 1 


90 00 


1 18 42 


32 


-4 00 


2 00 


-362±247 


-1.47 


nRR07 1 01 OR 


9Q37QS 


-1 70 


17 24 


64 


9 00 


9 00 


404 4- 1 SQ 


9 S4 


HRR07I 090 


294835 


-3 22 


1 14 


4 


-0 20 


40 


35±13 


2.69 


HRR0803 1 QR 
vJIVDVjoVjj 1 * 15 


306757 


-2 85 


S7 S7 


4 


-0 10 


14 


93-1-6 


3 83 

J.OJ 


nRR0803iQr" 

VJrvlJVJOVJ J 1 / V . 


306778 


-0 77 


1 3.3 1 


32 


-1 00 


1 00 


174±91 


1 .91 


HRR08041 1 


30001 o 


38 46 


48 4S 


4 


SO 
-u. JU 


SO 


1 1 6-I-9S 


4 64 


(TRR08041 3 A 

VJ 1\D VJ OVJ^ 1 J A 


309096 


-0.42 


9.05 


<^ 


-1 00 


1 00 


107±59 


1.81 


HRR08041 3R 

VJ 1\D VJ O yj'-r 1 JD 


3091 1 1 


-1 44 


4.96 


32 


-1.00 


1.00 


1 15±50 


2 30 


HRR080430 

VJ 1YU VJ O VJt 1 J 1/ 


310613 


-1.24 


12.84 


256 


-4.00 


4 00 


91 ±431 


21 


VJ IVU \J O \J VJ \J J U 


313087 


-0.54 


5.10 


16 


-1.00 


1.00 


5±59 


0.08 


nRR08060S 

VJIYU VJOV/VJVJ J 


313299 


-5.46 


15.53 


<^ 


-0 20 


20 


35±18 


1.94 


HRR080607 

VJIYU VJO VJVJVJ / 


313417 


-6.13 


12.05 


<^ 


-0 50 


50 


26±30 


0.87 


(TRR080791 

VJ 1\D \J O VJ 1 Zl 


317508 


-3 39 


8.64 


64 


-2 00 


2 00 


-86il 10 


-0.78 


HRR080Q1 f\ A 


3948QS 


-2 66 


39 58 


1 28 


9 00 


4 00 


<0^-L7 1 A 
jo J^Z, 11 


9 73 


frRRnSI 999 

VJIYUVJO 1 Z.Z.Z. 


337914 


-0 80 


15 58 


4 


- 1 .00 


1 00 


227±51 


4.45 


HRR0Q0494 
vjkijVJ :/^J i +i- i + 


3503 1 1 


-0 94 


4 95 


16 


90 


20 


14il4 


1 00 


^;RR0Q0^^1 8 

VJIYU VJ ~VJU 1 O 


3^083 


46 01 


135 35 


g 


-2.00 


2 00 


267 ±72 


3 71 


GRB090715B 


357512 


-4.80 


21.06 


16 


-2.00 


3.00 


275±155 


1.77 


GRB090812 


359711 


-6.93 


41.20 


256 


-6.00 


6.00 


-22±202 


-0.11 


GRB090926B 


370791 


-22.00 


36.00 


512 


-10.00 


8.00 


746±627 


1.19 


GRB091018 


373172 


-0.28 


2.92 


64 


-2.00 


1.00 


143±297 


0.48 


GRB091020 


373458 


-2.54 


13.84 


128 


-3.00 


2.00 


-187±177 


-1.06 


GRB091024 


373674 


-9.58 


27.29 


512 


-10.00 


10.00 


912±604 


1.51 


GRB091029 


374210 


-4.03 


38.98 


256 


-10.00 


10.00 


-112±395 


-0.28 


GRB091208B 


378559 


7.66 


10.61 


64 


-1.00 


1.00 


105±66 


1.59 


GRB 100621 A 


425151 


-6.79 


40.31 


256 


-3.00 


3.00 


1199±311 


3.86 


GRB100814A 


431605 


-4.40 


29.39 


256 


-4.00 


4.00 


862±147 


5.86 


GRB 100906 A 


433509 


-1.49 


26.16 


128 


-2.00 


2.00 


105±79 


1.33 


GRB 11 0205 A 


444643 


118.89 


293.99 


64 


-1.00 


1.00 


-29±52 


-0.56 


GRB110213A 


445414 


-3.42 


5.29 


512 


-3.00 


3.50 


602±746 


0.81 



Table 6. Correlation coefficients of the lag-iJpk relation. 



Coefficient Type 


Correlation Coefficient 


Null Probability 


Pearson's r 


-0.57±0.14 


4.83 X fCT 3 


Spearman's r s 


-0.50±0.12 


1.36 x lfj- 2 


Kendall's r 


-0.37±0.14 


1.18 x icr 2 



physical interpretation for spectral lag and a lag-luminosity rela- 
tion. In the literature, several possible interpretations have been dis- 
cussed (|Dermed[l 998; Salmonson 200(i| ;|loka & Nakamurall200ll; 
Kocevski & Liand |2003l: ISchaefej 12004; bin et al.l |2004 iRvdel 
2005l ; IShen et al.ll2005l ; lLu et al.ll2006l ; |Peng et alj201lh . 

One proposed explanation for the observed spectral l ag is the 
spectral evolution during the prompt phas e of the GRB jPermed 
1998; Kocevski & Lian j200ilRvdel2005l) . Due to cooling effects, 



Spk moves to a lower energy channel after some characteristic 
time. When the peak energy (i? p k) moves from a higher energy 
band to a lower energy band, the temporal peak of the light curve 
also moves from a higher energy band to a lower one, which results 
in the observed spectral lag. In a recent study, |Pengetal.ll201ll) 
suggest that spectral evolution can be invoked to explain both pos- 
itive and negative spectral lags. Hard-to-soft evolution of the spec- 
trum produces positive spectral lags while soft-to-hard evolution 
would lead to negative lags. In addition, these authors also suggest 
that soft-to-hard-to-soft evolution may produce negative lags. 

A schematic diagram showing a hard-to-soft scenario is de- 
picted in Fig. [5] Initially, Epk of the spectrum is in the high-energy 
band, which results in a pulse in the light curve of the high en- 
ergy band. Then E p ^ moves to the lower energy band resulting in 
a pulse in the low-energy light curve. The temporal difference be- 
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Table 5. Observer-frame and source-frame slopes and correlation coefficients of the lag-L; so relation. Conservative 10% 
uncertainty is assumed for cases without uncertainties. 



Energy Bands 


Frame 


Slope 


Correlation Coefficient 


JNumber or LrKBs 


Reference 


(0.3-1), (3-10) keV 


Observer 


0.95±0.23 




9 


Mareutti et al. (2010) 


(6-25), (50-400) keV 


Observer 


1.16±0.07 


-0 79+ 16 


8 


Arimoto et al. (2010) 


(15-25), (25-50) keV 


Observer 


1.4±0.1 


- 0.63±0.06 


21 


U10 


(15-25), (50-100) keV 


Observer 


1.5±0.1 


- 0.60±0.06 


28 


U10 


(15-25), (100-200) keV 


Observer 


1.8±0.1 


- 0.67±0.07 


27 


U10 


(25-50), (50-100) keV 


Observer 


1.2±0.1 


- 0.66±0.07 


27 


U10 


(25-50), (100-200) keV 


Observer 


1.4±0.1 


- 0.75±0.07 


25 


U10 


(25-50), (100-300) keV 


Observer 


1.14±0.1 




6 


Norrisetal. (2000) 


(25-50), (100-300) keV 


Observer 


0.62±0.04 


-0.72±0.07 


6 


Hakkila et al. (2008) 


(50-100), (100-200) keV 


Observer 


1.4±0.1 


- 0.77±0.08 


22 


U10 


(20-100), (100-500) keV 


Source 


1.23±0.07 


-0 90+ 12 


8 


Arimoto et al. (2010) 


(100-200), (300-400) keV 


Source 


0.9±0.1 


- 0.76±0.06 


22 


Ukwatta et al. (2010b) 


(100-150), (200-250) keV 


Source 


1.2±0.2 


- 0.82±0.05 


24 


This work 



1000 



> 

+ 

LU 



100 




T 

100 

Lag/(1+z) msec 



Figure 4. The source-frame peak energy (-E p k(l + z)) versus source-frame spectral lags. The energy bands, 100 — 150 keV and 200 — 250 keV, corresponding 
to the lag extractions are shown in hashed red bands on the plot. 



tween the two pulses in the light curves would then be a measure 
of the cooling time scale of the spectrum. 

If this were the only process that caused the lag then in a sim- 
ple picture one would expect the source-frame average U p k to lie 
within the two energy bands in question. According to Fig. [4] for 
the majority of bursts the source-frame E p k lies outside the energy 
band 100 — 250 keV, indicating that the simple spectral evolution 
scenario described above may not be the dominant process respon- 
sible for the observed lags. However, it is worth noting that a pulse 
in a specific energy band may not always mean that the E p ^ is also 
within that energy band. There are other issues associated with this 
model: 1) the calculated cooling times based on simple synchrotron 
models are, in general, relatively small compared to the extracted 
lags, and 2) short bursts which exhibit considerable spectral evolu- 
tion do not show significant lags. 



Another model that purports to explain spectral lags is based 
on the curvature effect, i.e., a kinematics effect due to the observer 
looking at increasingly off-axis annulus areas relat i ve to the line- 
of-sight Jsaimonson I2odo1;lloka & Nakamurall200ll ; lDermeil2004l 
Ishen et al .1120051 : iLu et al .1 120061) . Fig. [6] illustrates how the spec- 
tral lag could arise due to the curvature effect of the shocked shell. 
Due to a smaller Doppler factor and a path difference, the radia- 
tion from shell areas which are further off axis will be softer and 
therefore lead to a lag. As with spectral evolution models, there are 
difficulties associated with the curvature models too. These kine- 
matic models generally predict only positive lags. As can be seen 
from Table[3]some of the measured lags are negative, and therefore 
these lags present a real challenge for the simple curvature models. 

It i s possible that sp ectral lags are caused by multiple mech- 
anisms. IPeng et all d201ll) investigated spectral lags caused by in- 
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Time Evolving E pk 




Time 

Figure 5. The time evolution of the i? pk across energy bands may cause the observed spectral lags in GRBs. 



Hard 

Relativistically Expanding Shell \ 




Figure 6. Spectral lags could arise due to the curvature effect of the shocked shell. At the source, the relativistically expanding shell emits identical pulses from 
all latitudes. However, when the photons reach the detector, on-axis photons get boosted to higher energy (hard). Meanwhile off-axis photons get relatively 
smaller boost and travel longer to reach the detector. Thus these photons are softer and arrive later than the on-axis photons. 



trinsic spectral evolution and the curvature effect combined. They 
showed that the curvature effect always tends to increase the ob- 
served spectral lag in the positive direction. Even for cases with 
soft-to-hard spectral evolution, when the curvature effect is intro- 
duced lags become positive. Hence they predict that the majority of 
measured spectral lags should be positive, which is consistent with 
the findings of this work and U10. 



5 SUMMARY AND CONCLUSION 

We have investigated the spectral lag between 100 — 150 keV and 
200—250 keV energy bands at the GRB source-frame by projecting 
these bands to the observer-frame. This is a step forward in the 
investigation of lag-luminosity relations since most of the previous 
investigations used arbitrary observer-frame energy bands. 

Our analysis has produced an improved correlation between 



spectral lag (r) and isotropic luminosity over those previously re- 
ported with the following relation: 



r/(ms) 
1 + 2 



L iso /(ergs 1 ) 



10 54.7 



(8) 



We also find a modest correlation between the source-frame 
spectral lag and the peak energy of the burst, which is given by the 
relation, 



r/(i 



1 + z 



E pk (l + z)/(keV) 

1Q3.7 



(9) 



Finally, we mentioned two simple models and noted their lim- 
itations in explaining the observed spectral lags. 
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